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Abstract 

An attempt to study the feasibility of a new experiment to search for double beta decay in ^^^Sn 
and ^^^Sn was carried out by using ultra-low background HPGe detector (244 cm^) inside the 
Gran Sasso National Laboratory (LNGS) of the INFN (Italy). A small sample of natural Sn 
was examined for 2367.5 h. The radioactive contamination of the sample has been estimated. 
The data has also been considered to calculate the present sensitivity for the proposed search; 
half-life limits ^ 10 ^^ — 10 ^^ years for /3~^EC and EG-EC processes in ^^^Sn and ^ 10 ^^ years 
for transition in ^^^Sn were measured. In the last section of the paper the enhancement 

of the sensitivity for a proposed experiment with larger mass to reach theoretically estimated 
values of half-lives is discussed. 

Keywords: Double beta decay, ultra low background, HPGe 7 detector 


1. Introduction 

The existence of non-zero mass of neutrino has been established by neutrino flavor-oscillation 
experiments [1-3]. However, the nature of the neutrino, either Dirac or Majorana, can only 
be tested through observation of neutrino-less double beta decay, which is a rare second order 
transition involving two isobars. Moreover, the double beta decay experiments have the potential 
to establish the absolute scale of the neutrino mass, to prove the hierarchy of neutrino mass, to 
test the existence of right-handed admixtures in the weak interaction, and to test some other 
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effects beyond the Standard Model. The decay can occur as : 


(^4, Z) (A, Z 2') -\- 2e + (2z/e) (1) 

(^4, Z) [A, Z — 2) 2e^ + (2z/e) (2) 

e + (^4, Z) (^4, Z — 2) e"^ + {‘2.1^0) (3) 

2e + [A, Z) [A, Z — 2) (2z/e) (4) 


Triggered by the important implication of neutrino mass, a new generation of experiments 
are aimed to observe Oupp decay in various isotopes and with different experimental techniques 
[4]. Among the main experimental activities, we remind that recently three {3(3 experiments have 
published new experimental data: GERDA [5], EXO-200 [6] and KamLAND-Zen [7]. GERDA 
is searching for decay in ^^Ge while EXO-200 and KamLAND-Zen are looking for the 

decay in ^^^Xe. Several other experiments are going through their R&D phase. These include 
GUORE (i^OTe) [8], SuperNEMO (^^Se) [9], MAJORANA (^^Ge) [10], SNO+ (^^^Te) [11], NEXT 
(^^^Xe) [12] etc.. Present experiments are sensitive to half-lives of the order of 10^^ years and 
a corresponding effective mass of the neutrino ^ 100 meV [13]. Gonstrained by the 

uncertainties in the calculation of nuclear transition matrix elements [14], phase space factor [15] 
and by the qa value [16] that are used to determine {m^^) as well as by the different background 
for different isotopes, it is essential to measure the half-life for Oupp decay in several isotopes. 

Although Sn was considered as one of the potential candidates since the 1980s [17], there 
have been only a few studies on Sn. It can be mentioned that one of the earliest attempts for 
the experimental study of double beta decay was done by Eireman [18] in 1949 using ^^^Sn. 

Natural tin contains three isotopes which can decay via double beta transition; 
through two electron mode and ^^^Sn through ^d+EG and EG-EG processes, as given in the 
following: 


i22sn ^ i22Te + 2e“ + (2Pe) 

( 5 ) 

+ 2e“ + (29^) 

( 6 ) 

+ ^ + e+ + {2ve) 

( 7 ) 

2e-+ 

(8) 


The Q-values for the decay transitions in ^^^Sn, ^^^Sn and ^^^Sn are 1919.82T0.16 keV, 372.9T2.7 
keV and 2291.1±1.5 keV [19], respectively. The natural isotopic abundances for these isotopes 
are 0.97(1)%, 4.63(3)% and 5.79(5)% [20], respectively. The gamma rays from ^^^Sn decays are 
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at low energy where the background is high and makes it more difficult to study this decay modes 
in an external source experiment. Therefore, we have not considered this nucleus for the present 
study. 

Recently, investigations were carried out using natural tin samples for decay of ^^^Sn 

to the excited states of the daughter nucleus ^^^Te and and EC-EC processes in ^^^Sn by 

Dawson et al [21, 22], Kim et al [23] and Barabash et al [24]. Half-life limits of the order of 
10 ^^ — 10 ^^ years were obtained in those experiments. Searches for ^d+EC and EC-EC processes 
in ^^^Sn were also carried out by Kidd et al [25] and Barabash et al [26, 27] using enriched 
material obtaining the best half-life limit of 10^^ years. The Kims collaboration [28] gave a 
half-life limit 2.0 xlO^^ years for neutrino-less double beta decay of ^^^Sn using tin-loaded liquid 
scintillator. 

An attempt to study the feasibility of a new experiment to search for double beta decay in 
^^^Sn and ^^^Sn was carried out using an ultra-low background HPGe detector. The gamma 
rays produced by the de-excitations of the excited levels of either ^^^Cd or ^^^Te can be detected 
by the HPGe detector. The aim of the present study was to investigate the decay processes 
to the ground state as well as to the excited states in ^^^Sn and to the excited states in ^^^Sn 
using 7 -ray spectrometry. We will discuss in the last section also a possible enhancement of the 
sensitivity with a proposed experiment with larger mass to reach the half-life values estimated 
by theory. 

2. Experimental set-up and measurements 

A coaxial closed-end n-type ultra low background HPGe (GeBer) 7 detector (244 cm^) was 
used for the measurement. The detector is placed in the Gran Sasso National Laboratory (LNGS) 
of the INEN, Italy, located underground at ~ 3600 meters of water equivalent. The schematic 
diagram of the experimental setup is shown in Eig. 1. The physical dimensions of the detector 
are given in Table 1. More details about the detector and the ultra low background setup can 
be found in [29]. 

Monte-Carlo simulation of HPGe detector has been performed using the GEANT4 software 
library [30]. To validate the results of the simulation, we have compared the simulated efficiencies 
with known experimental results. The full-energy-peak (EEP) efficiency of the HPGe detector 
was measured for different 7 ray energies using ^^^Am, ^^^Ba, ^^^Cs and ^^Co point sources. 
The sources were placed at a distance of 26.5 cm above the end cap of the detector. The energy 
resolution (EWHM) of the detector was measured as 2.0 keV at 1332.5 keV. 
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Figure 1: Schematic diagram of experimental set-up used in simulation. Figure not in scale. 


A natural tin (purity 99.997%) sample of 13.3 g and with a thickness of 4.5 mm was placed 
on the end-cap, and the spectrum was measured with the HPGe detector in order to study the 
double beta decay processes in ^^^Sn and ^^^Sn nuclei. Data were accumulated for 2367.5 h with 
the sample and for 6109.4 h without the sample (background). The energy spectra of sample 
and background, normalized to the time of measurement of sample, are shown in Fig. 2. 

3. Results and Discussion 

3.1. Efficiency of the detector 

The FEP efficiency of the detector was simulated for the calibrated point sources using 
Monte-Carlo based software libraries GEANT4 [30] and EGS4 [31]. Generally, the results of 
the Monte Carlo simulations deviate significantly (> 10%) from the experimental results [32- 
38]. The difference can be due to two reasons: the uncertainties associated with the detector 
shape parameters provided by the manufacturer or incomplete charge collection in the crystal 
during the measurement process. Usually, all the information about the inner components of 
the germanium detector are not supplied by the manufacturer. Thus uncertainties arise due to 
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Figure 2: Energy spectrum with 13.3 g of natural tin sample (Sn) for 2367.5 h of measurement in comparison 
with background spectrum (Bg) of ultra low-background HPGe detector measured for 6109.4 h. The energy of 
the 7 lines are in keV. 


Table 1: The geometrical parameters of the detector. It is worth noting that the only difference between optimized 
and nominal values of the detector is the hole inner dead layer (0.3 mm nominal value). See text. 


Detector parameter 

Optimized value (MG2)(mm) 

Ge crystal radius 

30.85 

Ge crystal length 

81.70 

Garbon (G) window thickness 

0.76 

Ge crystal-C window distance 

4.00 

Hole radius 

4.95 

Hole length 

73.40 

Ge side dead layer 

0.0 

Ge front dead layer 

0.0 

Hole inner dead layer 

5.00 


possible additional absorbers or on the nature of the materials allocated around the Ge crystal. 
Furthermore, the dimensions provided by the manufacturer correspond to the time of assembly 
of the detection system at room temperature. Then changes in the mechanical support of the 
crystal due to contractions at low temperature may lead to changes in the detector configurations 
[33]. There is also a certain uncertainty in the crystal parameters e.g. the dead layer thicknesses 
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Table 2: Comparison of experimental efficiency with simulated efficiency computed with two different Monte- 
Carlo codes: EGS4 (E4) and GEANT4 (G4). Efficiencies related to MCI and MC2 are calculated considering 
two different thickness of the inner dead-layer: 1.5 mm and 5.0 mm, respectively. 


Source 

Energy 

Exp 

MCI 

Exp 

MCI 

MC2 

Exp 

MC2 

MCI 

Exp 

MCI 

MC2 

Exp 

MC2 


keV 

% 

(E4)% 

(E4) 

(E4) % 

(E4) 

(G4)% 

(G4) 

(G4)% 

(G4) 

241 Am 

26.3 

0.0048 

0.0064 

75% 

0.0064 

75% 

0.0066 

73% 

0.0064 

75% 


59.5 

0.0938 

0.1089 

86 % 

0.1088 

86 % 

0.1135 

83% 

0.0974 

96% 

i33Ba 

81.0+79.6 

0.0913 

0.1073 

85% 

0.1044 

87% 

0.1033 

88 % 

0.0979 

94% 


276.4 

0.0111 

0.0131 

85% 

0.0114 

97% 

0.0132 

84% 

0.0116 

96% 


302.9 

0.0270 

0.0305 

89% 

0.0271 

100 % 

0.0313 

86 % 

0.0280 

96% 


356.0 

0.0798 

0.0926 

86 % 

0.0799 

100 % 

0.0955 

84% 

0.0810 

99% 


383.8 

0.0110 

0.0126 

87% 

0.0107 

103% 

0.0125 

88 % 

0.0109 

101 % 


661.7 

0.0668 

0.0818 

82% 

0.0679 

98% 

0.0948 

71% 

0.0690 

97% 

® 0 Co 

1173.2 

0.0527 

0.0652 

81% 

0.0531 

99% 

0.0680 

78% 

0.0520 

101 % 


1332.5 

0.0478 

0.0600 

80% 

0.0476 

100 % 

0.0636 

75% 

0.0480 

100 % 


[34] or the distance between end cap and Ge crystal. It has also been observed, that there can 
be a substantial increase of the dead layer thickness after some years of operating time [35]. 
Therefore, to match the experimental efficiency values, the detector parameters and the dead 
layer thickness can be varied [38]. In the first step, the detector parameters provided by the 
manufacturer was used in the simulation. A fine adjustment was then made by varying only 
the inner dead layer thickness in a systematic way to match the experimental efficiency values. 
In order to show the dependence of the calculated efficiencies on the inner hole dead layer, the 
results obtained by considering two different thicknesses (1.5 mm and 5.0 mm) are compared with 
the experimental values in Table 2 and in Fig. 3. In the latter case the FEP efficiencies have 
been divided by the branching ratios of the corresponding 7 rays. From Table 2, we can see that 
for the optimized dead layer of the inner hole, thickness of 5.0 mm, the simulated efficiencies (for 
both EGS4 and GEANT4) are in very good agreement with the experimental results except for 
the very low energy region. The uncertainty is within 4%, which is quite reasonable considering 
the statistical fluctuations. However, at very low energy, difference is quite large. This could be 
due to the very small mean free path of 7 rays at this energy scale. Therefore, some 7 rays are 
absorbed before reaching the detector. The optimized parameters were then used to calculate 
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Figure 3: Comparison of the experimental efficiencies with those obtained from Monte-Carlo simulation using two 
different codes EGS4 and GEANT4. Here the FEP efficiencies have been divided by the branching ratios of the 
corresponding gammas. 


the FEP efficiency for a volume source of dimensions equal to the tin sample used in the present 
work. 

3.2. Radioactive contamination in the sample 

The sources of background radiations measured in an underground laboratory can be classified 
into several categories; environmental radioactivity including radon, gamma rays, neutrons from 
natural fission and from the ((a,n) reaction, radioactive impurities in the detector and shielding 
materials and cosmic rays with relevant contributions from muons and neutrons [39, 40]. In the 
measured spectra, background components can be identified by their characteristic 7 -emission 
peaks. The specific activities of the radioactive nuclei present in the natural tin sample were 
calculated with the formula [29] 

^ _ {^s/ts ~ Slj/tl)) 
mr] 

where Ss and Sb denote the area under a peak in the sample and background spectra, respec¬ 
tively. The measurement time of the sample and background spectra are denoted by tg and 4 , 
respectively. The mass of the sample is represented by m and rj is the efficiency of the full-energy 
peak detection (7 also takes into account the decay fraction of the 7 ray which was obtained 
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Table 3: The radioactive contamination present in the Tin sample measured with HPGe detector. Limits are 
given at 90 % C. L. 


Source 

Energy 

Count rate (mBq) 

Count rate (mBq) 

V 

Activity in 


(keV) 

from Signal 

from background 

(%) 

(mBq/Kg) 

228Ac 

911.2 

0.014T0.002 

0.012T0.001 

1.16 

< 30 

212pb 

238.6 

0.053T0.005 

0.048±0.003 

6.94 

< 15 

208rpj^ 

583.2 

0.021T0.002 

0.015T0.001 

5.83 

7±3 

212Bi 

727.3 

0.004T0.001 

0.006T0.001 

0.37 

< 34 

214Bi 

1764.5 

0.098T0.001 

0.006T0.001 

0.41 

42T22 

214pb 

351.9 

0.048T0.004 

0.033T0.002 

3.81 

30±9 

210pb 

46.5 

0.335T0.008 

0.196T0.004 

1.30 

801T51 

234Th 

63.3 

0.095T0.005 

0.103T0.004 

1.11 

< 28 

226 Ra 

186.2 

0.056T0.006 

0.060T0.004 

0.72 

< 79 

234mpg^ 

1001.0 

0.004T0.001 

0.003T0.001 

0.03 

< 739 


661.7 

0.030T0.002 

0.003T0.001 

5.03 

39T4 

6OC0 

1173.2 

0.020T0.002 

0.012T0.001 

3.63 

17±4 

40k 

1460.8 

0.094T0.003 

0.087T0.002 

0.34 

155T79 


from National Nuclear Data Center (NNDC) [41]). Activities of the natural radioactive sources 
present in the tin sample are listed in Table 3. Activity limits are given at 90 % confidence level 
(C. L.) according to the Feldman-Cousins method [42]. 


3.3. Double beta decay study 

In the accumulated spectra with the tin sample, no peaks are observed which could be 
unambiguously attributed to the double beta decay processes of Therefore, only half- 

life limits are calculated using the formula 


Ti/2 > 


(In 2 ) • • 77 • t 

lim S 


( 10 ) 


where N is the number of yd^d-active nuclei, t is the measurement time and 77 is the FEP detection 
efficiency. The expression lim S is the number of events of the effect searched for which can be 
excluded at a given C. L.. All the limits reported in the present study are given at 90% C. L.. 
The values of lim S were calculated using the Feldman-Cousins procedure [42]. The detection 
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efficiencies of the double beta processes in the tin isotopes were calculated using the GEANT4 
software library [30]. The DECAYO [43] event generator has been used to generate the initial 
kinematics of the particles. 

Taking into account the isotopic composition of natural tin, the sample contained 6.55 x 10^^ 
nuclei of ^^^Sn and 3.90 x 10^^ nuclei of ^^^Sn. The measured limits on half-lives for these two 
isotopes along with today’s best experimental limits and theoretical estimates are reported in 
Table 4. 
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Figure 4: Partial decay scheme of [44]. The energies of the excited levels and the emitted 7 quanta are in 

keV with relative intensities of 7 quanta are given in parenthesis. 


3.3.1. EC-EC decay of^^^Sn 

In the 2 z/ EC-EC decay of ^^^Sn to the ground state of ^^^Cd, all the energy release is carried 
away by the neutrinos except for a very small amount emitted as X-rays. These X-rays lie 
below the energy threshold of the measurement apparatus. In case of the neutrino-less mode, 
bremsstrahlung 7 quanta are emitted with an energy equal to E^ = Q 2 y — ci — €2 — E^^e^ where 
€i are electron binding energies of daughter nuclide, and Eexe is the populated level energy of 
^^^Cd. The partial decay scheme of ^^^Sn is shown in Eig. 4 [44]. Eor the transition to the 
excited state, the bremsstrahlung 7 quanta are accompanied by 7 rays emitted from nuclear 
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de-excitation. We did not observe any peak with the expected energies for the EC-EC decay of 
^^^Sn. Only lower limits of half-lives are obtained using the Eeldman-Cousins prescription [42]. 
The Ti /2 limits along with the 7 energies and the corresponding detection efficiencies are shown 
in Table 4. 




Figure 5: Part of the energy spectra accumulated with Sn sample for 2367.5 h of measurement by ultralow- 
background HPGe 7 detector for energy range (450-600) keV (upper panel). The energy spectrum accumulated 
without the sample over 6109.4 h is also shown (lower panel). Fits of the 511 keV annihilation 7 peaks are shown 
by solid line. 


3.3.2. EC processes in 

The {Ou + 2z/)/3+EC transition to the ground state of daughter nuclei is accompanied by 
two annihilation 7 quanta with energy 511.0 keV. Moreover, the detector surroundings and the 
sample contains the radioactive element ^^^Tl which emits 7 quanta of energy 510.8 keV. The 
energy spectra accumulated with and without the sample in the energy interval 450-600 keV 
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are presented in Fig. 5. The measured area of the annihilation peak is (130±7) counts. The 
area of the peak in the background spectrum is (100T15) counts (normalized on the time of 
measurement with the tin sample). The excess of events in the data accumulated with the tin 
sample can be explained by radioactive contamination of the sample itself. In particular, the 
decay of ^^^Tl contributes 17±7 counts of 510.8 keV gamma quanta. The difference in the areas of 
the annihilation peak (13T18) counts, which can be attributed to electron capture with positron 
emission in ^^^Sn, gives no indication on the effect. In accordance with the Feldman-Cousins 
procedure, 42.6 counts can be excluded at 90% C.L. Taking into account the FEP efficiency ( 77 ) 
of 7.58%, the lower limit of the (Oz/ + 2z/);d+EC transition to the ground state is 2.18x10^^ yr. 

In the case of the (0z/-|-2z/)/3+EC transition to the excited 2 ^ (617.5 keV) states of ^^^Cd, no 
peak has been observed. The limit on the half-life is reported in Table 4. The theoretical T 1/2 
is estimated to be very high because of the low phase space available for this transition. 

3.3.3. The (3~ decay of^‘^^Sn 

Eor decay of ^^^Sn, only transitions to excited states of ^^^Te can be studied as we are 

measuring gamma quanta. The partial decay scheme of ^^^Sn is shown in Eig. 6 [45]. No peak 
has been observed at 602.7 keV, as shown in Eig. 7. 

To derive the lim S value for the 602.7 keV 7- line, a model of the background was built from 
two Gaussian (one to describe the effect, and the second one to take into account the gamma 
peak with energy 609.3 keV of plus first-degree polynomial to describe the continuous 

background. The best fit (y^/n.d.f. = 50.8/51 = 0.99) achieved in the energy interval 586- 
628 keV gives the area of the peak searched for 2±12 counts (at the energy 602.7 keV), which 
conservatively provides lim S = 22 counts. Considering the transition to the 2 ^ (2182.4 keV) 
state (studied here for the first time), the calculated efficiency to detect 7 quanta with energy 
602.7 keV is 4.06%, thus T 1/2 > 1.36x10^^ yr. The half-lives limits for the transition of ^^^Sn 
to other excited states of ^^^Te are reported in the lower part of Table 4. 

The overall result is that due to the small mass of the tin sample ( 13.3 g), the obtained limits 
are rather poor compared to the existing experimental data. Therefore, it is certainly worthwhile 
to repeat the measurement with a larger sample, may be even enriched in the isotopes of interest, 
in order to obtain more stringent limits. 
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Table 4: The experimental limits and theoretical predictions for /3+EC/ECEC decay in ^^^Sn and /3 /3 decay 
in I24sn. 


Transition 

Mode 

Energy of 
7 rays 
(keV ) 

V 

(%) 

LimS 
90% C. L 

T®J^(yr) at 90% C. L. 
Present Previous 

work work [24, 27] 

[46, 47] 

-i. 








/3+EC; g.s. 

Ou T 

511.0 

7.58 

42.6 

2.18x10^'^ 

0.92x10“ 

3.8x1024 

/3+EC; 2+ 617.5 

Ou T 2z/' 

617.5 

5.95 

7.8 

9.40 xlO^'^ 

7.02x10“ 

2.3xl002 

ECEC; K^L^-,g.s. 

Ou 

1889.1 

2.16 

2.1 

12 .68x10^’’ 

8.15x10“ 

- 

ECEC; 2+ 617.5 

Ou 

617.5 

4.59 

7.8 

7.25x101’' 

9.40x10“ 

- 

ECEC; 0^ 1224.4 

Ou 

606.9 

3.32 

40.2 

1 .02x101’' 

12 .86x10“ 

- 



617.5 

3.22 

7.8 

5.09x101’' 

- 

- 

ECEC; 2j 1312.3 

Ou 

617.5 

2.35 

7.8 

3.71x101’’ 

8.89x10“ 

- 



694.7 

2.11 

35.9 

0.72x101’’ 

- 

- 



1312.3 

0.86 

18.0 

0.59x101’’ 

- 

- 

ECEC; Oj 1433.2 

Ou 

617.5 

2.85 

7.8 

4.50x101’’ 

6 .86x10“ 

- 



815.8 

1.01 

32.0 

0.39x101’’ 

- 

- 

ECEC; 2j 1468.8 

Ou 

617.5 

2.06 

7.8 

3.26x101’’ 

6.46x10“ 

- 



851.1 

1.59 

6.0 

3.27x101’’ 

- 

- 



1468.8 

1.05 

5.2 

2.49x101’’ 

- 

- 

ECEC; Oj 1871.14 

Oi^ 

617.5 

4.76 

7.8 

7.52x101’’ 

13.43x10“ 

- 



1253.4 

2.72 

24.4 

1.37x101’’ 

- 

- 

ECEC; 617.5 


617.5 

6.06 

7.8 

9.58x101’’ 

11.94x10“ 

4.9x1028 

ECEC; 0]^ 1224.4 

2 v 

606.9 

4.26 

40.2 

1.30x101’’ 

16.25x10“ 

7.4x1024 



617.5 

4.51 

7.8 

7.13x101’’ 

- 

- 

ECEC; 2j 1312.3 

2 v 

617.5 

3.26 

7.8 

5.15x101’’ 

11.24x10“ 

1.9xl002 



694.7 

2.96 

35.9 

1.01 xlOi’’ 

- 

- 



1312.3 

1.25 

18.0 

0.85x101’’ 

- 

- 

ECEC; Oj 1433.2 

2 v 

617.5 

3.90 

7.8 

6.16x101’’ 

8.64x10“ 

- 



815.8 

1.39 

32.0 

0.54x101’’ 

- 

- 

ECEC; 2+ 1468.8 

2 v 

617.5 

2.80 

7.8 

4.42x101’’ 

8.19x10“ 

6 .2xl00i 



851.1 

2.15 

6.0 

4.42x101’’ 

- 

- 



1468.8 

1.42 

5.2 

3.37x101’’ 

- 

- 

ECEC; 0;^ 1871.14 

2 v 

617.5 

4.76 

7.8 

7.52x101’’ 

13.43x10“ 

5.4x10®'* 



1253.4 

2.72 

24.4 

1.37x101’’ 

- 

- 

124S„ ^ 124Te 








2+ 602.7 

0 j'+2j/+0 !2 x° 

602.7 

6.05 

22.0 

2 .02x101® 

9.1x1020 

4.8x102o 

I3-/3--, 2+ 1325.5 
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Figure 6 : Partial decay scheme of [45]. The energies of the excited levels and of the 7 quanta are in keV 

with relative intensities of 7 quanta are given in parenthesis. 


4. Double beta decay study with enhanced sensitivity 

The sensitivity of a future experiment can be enhanced by either increasing the mass of the 
sample and the efficiency of the detector or by reducing the background. The efficiency of the 
detector depends on the position and the shape of the sample. The 7 rays can undergo either self¬ 
attenuation to loose in the sample itself or absorption in other interposed media. The intensity 
of the gamma ray decreases to half of its initial value passing through a thickness of a half-value 
layer (HVL). The HVL depends on the atomic number of the element and also on the energy 
of the 7 ray. HVL for tin material at different 7 -ray energies were calculated from the mass 
attenuation coefficients (/i/p) values at different energies taken from NIST data [48]. 

Efficiencies of the detector were simulated for various sample thicknesses and different sample- 
detector distances. At a fixed 7 ray energy, though the efficiency initially increases for lower 
energy upto 400 keV and then decreases with an increasing thickness of the sample as shown 
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Figure 7: Part of the energy spectra accumulated with Sn sample for 2367.5 h of measurement. In particular, 
the region around the 602.7 keV peak, expected for the processes in to the excited states of 

is shown. Fits of the excluded peak at energy 602.7 keV along with peak due to with energy 609.3 keV are 

depicted. Peak of 208 energy 583.2 keV is also shown. 


in Fig. 8, the number of double beta nuclei increases with increasing the mass of the sample. 
Hence, the product of the efficiency and the number of nuclei together increases with increase in 
the mass of the sample as shown in Fig. 9. 

For a cylindrical sample of fixed radius, the growth of the FEP count tends to saturate at a 
certain sample thickness as shown in Fig. 9. With further increase in sample thickness, the FEP 
count alters very less but the Compton continuum due to high energy photons increases. As a 
result of the increased Compton background, the important low-intensity gamma-lines emitted 
due to double-beta decay processes may not be detected experimentally. Therefore, it is very 
important to optimized the sample thickness. The optimization of sample thickness for cylindrical 
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Figure 8: Detection efficiency of the HPGe detector at different energies in the region-of-interest for different 
thicknesses of the sample are shown. Legends are the thicknesses of the sample. 




Figure 9: Comparison of efficiency times mass of Figure 10: Correlation between sample thickness 
the sample for different thicknesses of the sample. and first derivative of efficiency times mass of the 

sample. 


geometry that maximizes the detection efficiency are studied by Barrera et ah [49], Shweikani et 
ah [50] and Li et ah [51] for HPGe detectors. For a cylindrical sample volume with fixed radius, 
optimal thickness has been determined by calculating the minimum value of the thickness that 
makes the first derivative of the sample mass multiplied by the detection efficiency nearly zero 
[50]. The correlation between sample thickness and efficiency times mass of the sample can be 
expressed by the Box Lucas function having the form [50] : 

7] ■ M = a{l - (11) 

where a, b are two constants and h is sample thickness. The variation of {rj ■ M) with sample 
thickness is shown in Fig. 9. In order to estimate the optimum sample thickness at a given energy, 
the correlation between the first derivative of efficiency times mass of the sample, (77 • M)' and 
sample thickness was studied. From Fig. 10, we can see that the value of (77 • M)' tends nearly 
to zero at sample thickness 36 mm and higher for gamma energies 600-2000 keV. We consider 
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= 36 mm as optimal thickness for the future large-scale experiment. 
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Figure 11: Schematic diagram of the experimental set-up of the proposed future experiment. Figure is not in 
scale. 

Considering the geometry of the present HPGe detector, a cylindrical shaped (086.0 mm x 
36.0 mm) tin sample of mass ^1530 g (density of tin is taken as 7.31 g/cm^) can be placed on 
the sample holder, as shown in Fig. 11. A tin sample, enriched in ^^^Sn up-to 90%, will contain 
7.40x10^^ nuclei of this isotope. Enriching in ^^^Sn instead by the same factor, the sample will 
contain 6.68x10^^ nuclei of it. Considering the efficiencies 2.44 % (2z/ EC-EC) and 2.36 % (Oz/ 
EC-EC) for decay of ^^^Sn to the 2^ (617.5 keV) state of ^^^Cd, the following sensitivity limits 
(at 90% C. L.) could be reached after two years of measurements: 

g.s. ^ 617.5 keV) > 3.0 x IO^Ve 

g.s. ^ 617.5 keV) > 3.0 x lO^Vr. 

This result would be better than the previously measured best half-life limit[27]. Eor the double¬ 
beta decay transition of ^^^Sn to the excited 2^ (602.7 keV) state of ^^^Te, with 90% enriched 
^^^Sn and an efficiency of 2.50 %, the sensitivity would be the same as for ^^^Sn, up to 10^^ yr 
(at 90% C. L.) for the half-life. 
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5. Conclusion 


In our preliminary study for the search of excited state transitions concerning possible double 
beta decays of tin isotopes with a HPGe detector, we measured half-life limits of 10^^ — 10^^ 
yr for ;d+EC and EC-EC processes in ^^^Sn and 10^^ yr for /3~ transition in ^^^Sn. Eor the 
Oz/EC-EC decay of ^^^Sn to the ground state of its daughter nuclide we obtained a half-life limit 
of 1.27x10^^ yr. We showed that an experiment with larger mass, and material enriched to high 
levels of either ^^^Sn or ^^^Sn, could reach with the same HPGe detector half-life limits for both 
tin isotopes on the order of 10^^ yr for the decays to the excited levels. In India, an effort has 
been started to build up a bolometric Sn detector for experimental study of 0uf3f3 decay in ^^^Sn 
[52]. The experimental setup will be housed at upcoming India-based Neutrino Observatory site. 
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